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The path integral by which quantum field theories are defined is a particular solution of a set 
of functional differential equations arising from the Schwinger action principle. In fact these 
equations have a multitude of additional solutions which are described by integrals over a 
complexified path. We discuss properties of the additional solutions which, although generally 
disregarded, may be physical with known examples including spontaneous symmetry breaking 
and theta vacua. We show that a consideration of the full set of solutions yields a description 
of phase transitions in quantum field theories which complements the usual description in 
terms of the accumulation of Lee- Yang zeroes. In particular we argue that non-analyticity 
due to the accumulation of Lee- Yang zeros is related to Stokes phenomena and the collapse of 
the solution set in various limits including but not restricted to, the thermodynamic limit. A 
precise demonstration of this relation is given in terms of a zero dimensional model. Finally, 
for zero dimensional polynomial actions, we prove that Borel resummation of perturbative 
expansions, with several choices of singularity avoiding contours in the complex Borel plane, 
yield inequivalent solutions of the action principle equations. 
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1 Introduction 



Many of the ideas presented in this paper have also appeared previously in [1]. The intent of this 
article is to assemble what we view as the important points in a short and coherent summary, and 
to add some results concerning a relation between Lee- Yang zeroes and Stokes phenomena. 

Quantum field theories may be defined either by a path integral or by a set of functional 
differential equations which follow from the Schwinger action principle [2], 

S{tM = {h\5S\t2) . (1) 

For the sake of illustration, consider a zero dimensional "quantum field theory" defined by the 
action S{(f)). The generating functional for correlation functions of is 

/ + 00 
rf0exp(-5(0) + J0) ; (2) 
-oo 

where it is assumed that the integral is convergent. This is a solution of the action principle 
equations 



where Qi are the parameters of the theory. For the specific action 



{S'idj) + J) Z[J) = (3) 
d,^-^)zig,J)=0; (4) 



5(0) = l/i02 + |0^ (5) 



equations (3) and (4) become 



{gdj + fLdj + J)Z{J) = (6) 
{d,-\dj)Z{J) = {d,-ldj)Z{J)=0. (7) 

Although all these equations are included in the Schwinger action principle, it will be convenient 
to refer to equations involving variations of the fields (e.g., (3)) as Schwinger-Dyson equations, and 
those arising from variation of parameters (e.g., (4)) as action principle equations. 

In general the equations (3) and (4) have a several parameter class of solutions. For the action 
(5), the corresponding equations (6) and (7) have a three parameter class of solutions, which 
includes (2). The Schwinger-Dyson equations are solved by 

Z{J) = y^cj{g,fi) [ #exp(-5(0) + J0); (8) 

where Tj are inequivalent integration paths in the complex plane over which the integral converges 
and c/ are arbitrary functions of the coupling constants g and fi. The number of such contours 
matches the order of the differential equation (3). Figure 1 shows the domains of convergence, 
cos(4arg(0)) > 0, and a basis set of contours for real positive g. 
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Figure 1. 



Roughly speaking, the action principle requires the coefficients c/ to be independent of the 
parameters fi and g. This statement is imprecise since, if g is is taken to be complex and the 
argument of g is varied sufficiently, the contours of integration Fj must be changed to maintain 
convergence. For a general polynomial action, this statement holds for the "top coupling" asso- 
ciated with the highest power of in the action. For a given top coupling g, each contour Tj 
belongs to an equivalence class of contours for which integration gives the same result. The equiv- 
alence classes (see figure 2) consist of contours for which |0| — >■ oo within the same domains of 
convergence; the action is analytic in except at infinity. There is always a choice of contour Tj 
within an equivalence class which can be held fixed while making infinitesimal variations of the 
top coupling. The action principle requires that the coefficients associated with these contours do 
not vary as one makes infinitesimal changes in the couplings. 



Figure 2. For arg(5f4) = 0, (a) the contours B and B' are equivalent. However for 
arg(5f4) = 2 7r/3, (b) the integral over B' remains convergent while the integral over B diverges. 

To dispel doubt that the "exotic" solutions with complex integration contours have physical 
relevance, consider the effective action for the theory defined by (5), with real < and real 
g > 0. The expectation value for is a solution of the equation 




a) 



B' 





J 



d4> 



(9) 



2 



where F is the one particle irreducible effective potential. The tree level effective potential (in this 
case the action) has minima at = 0, ^^/Jifg. The minimum (0) = \fjifg + ■ ■ ■ corresponds 
to a solution of the Schwinger-Dyson equations with an integral representation involving a sum of 
the contours A and C drawn in figure 1. These integrals have contributions only from the saddle 
point at = \/ —1^1 g- Hence symmetry breaking vacua, whose existence is normally attributed 
to a thermodynamic limit, exist even in zero dimensions when the full set of solutions of the 
Schwinger-Dyson and action principle equations are considered. Note that no symmetry breaking 
term has been added to the action; the symmetry is broken by the choice of integration path^. 

Note that the solutions associated with the contours A and C in figure 1. have complex parts 
which do not appear at any order in a perturbative expansion about the saddle point. However, due 
to the linearity of the Schwinger-Dyson equations, one can sum the contours to get a solution for 
which the non-perturbative contribution is real. The reader might be concerned that the "exotic" 
solutions which are not integrals over real do not satisfy the axioms of Euclidean quantum 
field theory. Although one can easily arrange for the Greens functions to be real, one might still 
worry that even Greens functions are not manifestly positive. In fact one should postpone these 
questions for the higher dimensional case as perversities of the exotic solutions could vanish in the 
thermodynamic/continuum limit. In fact, we expect this to be the symmetry breaking 

vacua and theta vacua are examples of exotic solutions. 

In the subsequent section, we argue that the appearance of phase boundaries in quantum 
field theories is intimately related to a collapse of the solution set in the thermodynamic limit. 
This proposal can be made concrete in a zero dimensional analogue, for which we demonstrate 
a correspondence between the collapse of the solution set as a top coupling is set to zero and 
the accumulation of Lee- Yang zeros, leading to a non-analyticity in a coupling constant. In this 
context, the limit of a vanishing top coupling is analogous to the thermodynamic limit. The 
complementary descriptions of phase boundaries in terms of the accumulation of Lee- Yang zeroes 
and the collapse of the solution set share a common origin in Stokes phenomena. 

In section 3, we prove the equivalence of Borel resummation of the perturbative expansion 
about saddle points and exotic solutions of the Schwinger-Dyson equations, for various singularity 
avoiding contours in the Borel plane. Finally, in section 4, we argue that the action principle may 
emerge from the Schwinger-Dyson equations under suitable conditions, due to the collapse of the 
solution set in the thermodynamic limit, rather than being an independent set of equations. 

2 Collapse of the Solution Set and the Accumulation of 
Lee- Yang Zeros 

When complex values of the parameters of a quantum field theory are considered, phase boundaries 
which appear in the thermodynamic limit can be understood in terms of the accumulation of 
zeroes of the partition function which pinch the real axis, known as Lee- Yang zeroes [3]. As 
we shall explicitly demonstrate below, the accumulation of Lee- Yang zeroes is not confined to the 
thermodynamic limit. In zero dimensional polynomial theories, Lee- Yang zeroes also accumulate in 
the limit that the top coupling goes to zero. In this limit, a non-analyticity in the coupling gT-\ 
appears. In this context, we will propose an alternative (and complementary) description of the 
appearance of phase boundaries in terms of the collapse of the solution set of the Schwinger-Dyson 

^The standard way to non-pcrturbatively describe symmetry breaking vacua is to introduce a small symmetry 
breaking term which is taken to zero only after taking an infinite volume limit. 
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and action principle equations. 

Consider the zero dimensional action S = YlJ=i9i^^- There is a branch point at gr = 0, due 
to the necessity of rotating the contours of integration in the integral representation, in order to 
maintain convergence as the argument of qt is varied. It is always possible to find a contour in 
an equivalence class which can be held fixed under infinitesimal variations of the top coupling (see 
figure 2), so as to satisfy the action principle. However large variations in the argument of the 
top coupling require changes in the contour. In particular, under a rotation by 2 tt T the solutions 
transits among T Riemann sheets. The solutions are analytic in the couplings gi for / < T, except 
at infinity, since these couplings do not effect the domains of convergence. In the limit qt — > 0, the 
solution develops a branch point in the new top coupling, gr-i- The hmit — ^ is analogous to 
a thermodynamic limit. For the case T = 3 we will explicitly show that this limit is accompanied 
by the accumulation of Lee- Yang zeroes in the complex Qt-i plane. 

The appearance of this branch point can be understood in terms of the collapse of the solution 
set of the Schwinger-Dyson and action principle equations. In the limit qt — 0, with fixed aTg{gT) 
solutions of the Schwinger-Dyson and action principle equations either diverge, vanish or have a 
finite limit. It is easy to see which by considering the overlap of the domains of convergence in the 
complex plane for the case in which qt is the top coupling or for which Qt-i is the top coupling 
{St = 0). If a convergent contour for ^ is equivalent to one which lies within a single wedge 
of convergence for the case qt = 0, then the integral will vanish in the Qt ^ limit with fixed 
aTg{gT)- On the other hand if a convergent contour for 7^ is not equivalent to any contour 
lying within the wedges of convergence for qt = 0, then the Qt ^ limit is divergent. A finite 
(?r ^ limit exists only if an equivalent contour lies within two different wedges of convergence 
for the Qt = case. Figure 3 illustrates the various possibilities for the case T = 4. 




Figure 3. For arg(5f4) = and arg(5f3) = 3 7r/4, the integral over the contour or the equivalent 
contour F^ vanishes in the ^ limit, while the integral over F^ is finite. For arglg^) = and 

arg((yf3) = TT the integral over Tq diverges. 

Suppose that the argument of qt is kept fixed a.s qt ^ 0- In this limit the behavior of the 
partition function, defined by a particular contour of integration in the complex (j) plane, will change 
discontinuously as the argument of Qt-i is varied across certain critical values. For example the 
Qt limit may go from finite to divergent upon crossing a Stokes line. However it is possible to 
keep this limit finite by considering variations in the contour of integration which violate the action 
principle by terms which vanish in the qt ^ limit. For qt = 0, these variations are equivalent 
to analytic continuation in Qt-i- 
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To illustrate how this works, consider the case T = 4. As one varies the argument of keeping 
that of (74 fixed, a generic solution of the Schwinger- Dyson and action principle equations will 
become alternately divergent, finite, or vanishing in the (74 — limit. One can can keep the limit 
finite by adding contours such as T a in figure 3 when the argument of enters a wedge in which 
the contribution from this contour vanishes as (74 — >■ 0. 

This process is illustrated in figure 4. The domains of convergence for (74 7^ are indicated 
in light gray, while the domains of convergence for (74 = 0, (73 7^ are indicated in dark gray. 
Starting with a solution of the Schwinger-Dyson and action principle equations corresponding the 
contour A, aig{g^) is varied from tt to 7r/2 keeping arg{g/^ = 0. Initially the solution is finite as 
g4 — »• 0, since the contour A lies asymptotically within two domains of convergence for (74 = 0. At 
arg{g3) = 3 7r/4 the contour is modified to C = A + B. Note that the integration over B vanishes 
in the (74 limit for a neighborhood of aiglg^) = 3 7r/4, since B is equivalent to a contour lying 
within a single domain of convergence when g^ = 0. One can continue varying aiglgs) to 7r/2 
without a change in the asymptotic behavior as (74 — 0; the partition function remains finite in 
this limit. On the other hand, had the contour been fixed as A, the integral would be divergent 
as 5^4 —s> for a.Tg{gs) on the other side of the Stokes line &Tg{gs) = Svr/S. For arg^g^) = 7i/2 and 
(74 = 0, the contour A does not lie within the domains of convergence asymptotically, unlike the 
the contour C. 




Figure 4. As the argument of g^ is varied, the contour is changed to keep the integral finite in the 

(74 — > limit. 

Considering larger variations of arg((73) one can repeat this process to give a solution of the 
Schwinger-Dyson equations which is finite and satisfies the action principle in the (74 — > limit. For 
g4 = this process corresponds to analytic continuation in g^. The changes in contour necessary 
to maintain a finite (74 ^ limit (and avoid Stokes phenomena) give rise to the third order branch 
point at (73 = 0. 

We thus arrive at the picture that non-analyticity in a coupling constant arises due to the 
collapse of the solution set of the Schwinger-Dyson and action principle equations. In the example 
we have given, non-analyticity in gr-i arises as the top coupling g-r —>■ 0. In this limit, some 
solutions are finite, while others diverge or vanish. Because the equations are linear and certain 
solutions vanish in this limit, various classes of solutions with inequivalent integral representations 
also coalesce as gr — ^ 0. This permits larger variations of the contour which satisfy the action 
principle. At the same time larger variations of the contour may become necessary to obtain a 
finite partition function as coupling constants are varied. While we have explicitly demonstrated 
this mechanism for the g^ ^ limit of a zero dimensional polynomial theory, we propose that 
the same phenomena hold true for the thermodynamic limit, N 00 where is the number 
of degrees of freedom, in a mult i- dimensional field theory. In other words, non-analyticity in a 
parameter of a quantum field theory should be related to the collapse of the solution set of the 
Schwinger-Dyson and action principle equations in the N ^ 00 limit. 
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The analogy we have proposed between a non-analyticity arising in the gr ^ hmit of a zero 
dimensional theory and non-analyticities arising in the thermodynamic limit is strengthened by 
noting that the gr limit is accompanied by the accumulation of Lee- Yang zeroes along Stokes 
lines in the complex Qt-i plane. One can explicitly see how this occurs for T = 3. Consider the 
partition function 

Z = [ rf0e-(f^'+^^'); (10) 
Jc 

for positive real g and a contour C which goes to infinity with arg(</)) = ±2 7r/3, as in figure 5. 
The integral (10) can be evaluated in terms of an Airy function; 

Z = 27rze-^#r^/Mz(^). (11) 



Figure 5. 

The zeroes of Ai{z) lie along the negative real z axis. Since the argument of the Airy function 
in (11) is 2; = /i^/4(7^/^, zeroes of the partition function accumulate on the imaginary axis in the 
complex fi plane as (7 — > 0, pinching the real axis at /i = 0. In fact, arg(/i) = ±7r/2 are Stokes 
lines. The g ^ behavior of Z in the neighborhood of arg{f^i) = — 7r/2 can be seen by inspecting 
figure 6, in which the domains of convergence for g ^ and g = are indicated by the lighter and 
darker shaded regions respectively. For arg(yu) = — 7r/2 + e, Z diverges in the g ^ limit, with 
the leading term in a saddle point expansion given by 

Z~./?^exp(-i4). (12) 
V /i 6 g^ 

On the other side of the Stokes line, arg(/i) = — 7r/2 — e, Z converges in the g — > limit, with 
the leading term in a saddle point expansion given by 
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Figure 6. On either side of the Stokes hne at arg(/i) = —'njl^ the integral is either finite (c) or 
divergent (a) in the — > hmit. Zeroes of the partition function accumulate as (7 ^ for values 

of /i on the Stokes line (b). 

On the Stokes line, arg(yu) = — 7r/2, the two saddle point expansions (12) and (13) become 
comparable in magnitude as (7 ^ 0, since the real part of the exponential in (12) vanishes. The 
integral over the contour in figure 5 is equivalent to the sum of the integrals over two constant 
phase (steepest descent) contours, C\ and in figure 7, which pass through classical solutions for 
which the real part of the action is degenerate. The accumulation of zeroes on the Stokes line is 
related to the fact that the relative phase of the two integrals oscillates wildly with variations of 
\i in the (7 — limit, due to the factor of exY>{l^^ 1 9^) = exp(i l/xp/^f^). 



> 



C'2 

Figure 7. Steepest descent contours passing through the classical solutions at = 0, —fJ^/g, on 

the Stokes line arg(yu) = — 7r/2, a.Tg{g) = 0. 

Thus, the accumulation of Lee- Yang zeroes and the collapse of the solution set of the Schwinger- 
Dyson and action principle equations are complementary descriptions of the appearance of non- 
analyticity in a parameter of a zero dimensional theory as the top coupling qx — > 0. The two 
descriptions have a common origin in Stokes phenomena. We conjecture that these are also com- 
plementary descriptions of non-analyticity in parameters of quantum field theories which arise in 
the thermodynamic limit N — »• 00. Upon completion of this work, we discovered [1], in which 
Lee- Yang zeroes in a ID model with a wetting transition were shown to lie along Stokes lines 
associated with the asymptotic expansion in A^. The correspondence between Stokes lines and 
Lee- Yang zeroes has also been suggested in [5], in the context of Ising and gauge models. 

3 Borel Resummation and Complexified Path Integrals 

It is well known that the loop expansion in quantum field theory is an asymptotic rather than a 
convergent series. One approach to obtain non-perturbative information from the loop expansion is 
the Borel resummation, in which a convergent series (the Borel transform) in a variable t is obtained 
from the asymptotic series in h. The Borel transform is then inverted to give a function having 
the correct asymptotic expansion, but which contains non-perturbative information. Since there 
are actually an infinite number of such functions, differing by essential singularities at ^ — >■ 0, it is 
a non-trivial statement that the Borel transformation corresponds to the path integral. Moreover, 
there are frequently singularities of the Borel transform on the positive real t axis, due to instantons 
and renormalons, which prevents an unambiguous Borel resummation. 
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Under generic conditions, the number of classical solutions is equal to the number of indepen- 
dent solutions of the Schwinger-Dyson equations. For an arbitrary polynomial action in a zero 
dimensions, 



T 



n=l 



n 



(14) 



we will prove that various Borel resummations of perturbative expansions about classical solutions 
satisfy both the Schwinger-Dyson equations and the action principle equations, and therefore 
correspond to various complexified path integrals. 
Consider the partition function 



>i5W ■ 



(15) 



where the path F is equivalent to a steepest descent path passing through a dominant saddle point 
(classical solution) </> = 0^. The loop expansion yields a contribution to the generating function of 
the form; 




(16) 



n=0 



This series is asymptotic, but its Borel transform defined by 



TX 



5"(0a)^F(n + i; 



(17) 



converges to 



Bo.{t) 



1 



2m fc '^t-(^(0)-^(0„)) ' ' ' 

where in the vicinity of t = the contour C encloses, in the opposite sense (see figure 8), the two 
poles (j)a,j{t), for j = 1,2, which coalesce to 0q, as t — * 0. All the other poles are taken to lie outside 
the contour. The Borel transform has a singularity when one of the exterior poles coalesces with 
one of the interior poles, which occurs when t = S{(j)a') — S{(j)a) where (pa' is a neighboring classical 
solution. Doing the integral gives 



(19) 



0O 



Figure 8. 
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Thus far everything we have said is well known [6]. We now exhibit an exact relation between the 
Borel resummation and the exotic solutions of the Schwinger-Dyson and action principle equations. 
We invert the Borel transform by writing 

Z, = e-^^(^") [ dte-T.^ (20) 

with an as yet unspecified integration contour S in the complex t plane. Equivalence of (20) with 
a solution of the Schwinger-Dyson equations requires 

£ ,,__A_^^ . * £ M (21) 

where F/ indicate open paths in the complex plane which asymptotically lie within the domains 
of convergence determined by the the top coupling qt- Since C and cjVj are inequivalent, (21) 
involves a non-trivial exchange in the order of integration under which contours are not preserved. 
Instead of directly proving (21), we will show that (20) solves the Schwinger-Dyson equations and 
satisfies the action principle. 

We set ^ = 1 in what follows. If Za satisfies both the Schwinger-Dyson equations and the 
action principle then it is annihilated by the operators; 

T g 

L = J2(n-l)g^- ; (22) 

n=2 9n-l 



and 



d 1 (9" 

^- = ir 1^- (23) 

ogn n dgl 

It is convenient to define the quantity 

To show that LZa = 0, for Za defined in (20), it suffices to show that = where 



t ^ 5^(n - l)gnD,„_, -g,- ^ — - — 5(0„) . (25) 

71=2 ^" 

Before proceeding, we list several simple but useful identities. Due to the equation of motion, 
S'icbn) = 0, one has 



^ -S{^a) = (26) 



dgn n 

Identities obtained by differentiating the relation t — S{(f)a,j{t)) + S{(j)a) = 0, which defines 

')a,j{t), are 

d , 1 



(27) 



9 ^ _ KC, 
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where the equations of motion have been used again in deriving the last equation. 
Using these identities and the equations of motion, one can show that 



- l)gnD,^_, - 9i]^i7—-^ = ; (29) 

n=2 ^'^"'J'^ 

implying = 0, or LZ^ = 0, for any integration path in t. 

We are not done yet, since the equation = is a combination of the Schwinger-Dyson and 
the action principle equations. The integration path E will be constrained further by requiring the 
action principle to be separately satisfied. To this end, consider the quantity 



A^[klD,,D,, + {k + l)D,,jF 



(30) 



which will vanish if the action principle is also satisfied. Using the same identities discussed above, 
A may rewritten as 



A 



E dt 



-iVkl 



d_d_ 

dgk dgi 



4>a,j{t) 



E(-i) 



'kl 



d_d_ 

dgk dgi 



(31) 



ar 

= +00, avoiding 

a,l{t) and (f)a,2(t) 



Clearly A vanishes if the contour S begins at t = and ends at Re{t) 
singularities. The contribution from the boundary at t = vanishes because 
coalesce at t = 0, so that the factor of ^j(— 1)-' in (31) leads to a cancellation. Closed contours 
encircling branch cuts in the complex t plane also give 71 = 0. So long as a contour E for which 
yi = is used, the Borel resummation gives a solution of the Schwinger-Dyson and action principle 
equations. 

It would be very interesting to see to what extent this analysis extends to theories with non-zero 
dimension. The analysis is surely more difficult since, among other possible complications, there 
are singularities in the Borel plane due to renormalons as well as finite action classical solutions 
(instantons). 

The "exotic" solutions of the Schwinger-Dyson equations given by Z = cjZi, where the Zj 
are generated from classical solutions by Borel resummation, may in some sense be thought of 
as a generalized form of theta vacua, in which the cj play the role of theta parameters. In the 
usual approach a particular theta vacuum is selected by adjusting a surface term in the action and 
integrating over real fields. The surface term term effects the Schwinger-Dyson equations only at 
the space-time boundary, so in the infinite volume limit its role is only to set a boundary condition. 
It does this by putting a different weight on the contributions to Z coming from perturbative 
expansions about different classical solutions. We conjecture that when appropriately resummed, 
this is equivalent to a weighted sums over complexified path integrals. 



4 Emergence of Action Principle in the Thermodynamic 
Limit 

So far we have treated the Schwinger-Dyson and action principle equations as independent. This 
is certainly true for a finite number of degrees of freedom. If the action principle is not imposed, 
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the manner in which a solution of the Schwinger-Dyson equations changes as one varies a coupling 
is un- determined, as the there is a continuous set of solutions to the Schwinger-Dyson equations 
for any value of the coupling. However, with certain assumptions, the action principle arises due 
to the collapse of the solution set of the Schwinger-Dyson equations in the thermodynamic limit. 

Suppose that the solution set collapses in the thermodynamic limit, i.e. some solutions coalesce 
while others diverge, such that are solutions of the form Z ~ exp(— A^?") as — oo where there 
are only discrete possibilities for the free energy J . Discreteness of the solutions automatically fixes 
the dependence on the couplings. Let us assume that this dependence is described by a differential 
equation of the form 

O^Z = {^- ka) Z[g^, J^] = 0; (32) 

where K is a linear operator^. We will show below that these equations are equivalent to the action 
principle equations. The argument follows from the commutation relations of operators associated 
with the action principle and the Schwinger-Dyson equations. 

Writing the action in the form S{(j)i} = gafa{4>i}, the operators associated with the Schwinger- 
Dyson equations are 



dfa 

goL 



- J^ (33) 

while those associated with the action principle are 

Ha = dg^ — /ali^^j^jg^^} . (34) 

Thus one obtains the commutation relations 

[L,,^,]=0, [H^,Hfs] = 0. (35) 
Let us now write Oa = + A^. If LiZ = and OaZ = 0, then [Oa, Li]Z = 0, or 

LiAaZ = . (36) 

If A^Z is non-zero, it is a solution of the Schwinger Dyson equations. Moreover there is 
only one discrete possibility: AaZ = Ca{g/3}Z, where Cq, is some function of the couplings, so 
that HaZ = —Ca{g/3}Z. The coefficients Cq, can be absorbed by a coupling constant dependent 
re-scaling Z ^ Z' = e^^^i^^Z, where 

HaZ' = [Ha, e^]Z - e^'caZ = ((9,„ - cje^) Z = 0. (37) 

The existence of a solution to the equations {dg^ — Ca)e^ = requires (9g^c^ — dg^Ca = 0, which 
follows from {dg^cp — dgi^Ca)Z = —[Ha,Hp\Z = 0. The re-scaled partition function satisfies both 
the Schwinger-Dyson and Schwinger action principle equations, L^Z' = HaZ' = 0, even though we 
started with just the Schwinger-Dyson equations. 



^Thc operator Ka is necessarily linear for this equation to make sense in the thermodynamic limit. 
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5 Non- Polynomial Actions 



Although we have focused on polynomial actions, our discussion of the solution set of the Schwinger 
action principle equations can be readily extended to non-polynomial actions. A simple non- 
polynomial action is that one plaquette QED, with action S = [3 cos 9. The generating functional 



Z{J, J) = 

is a solution of the differential equations 



d9e 



/3cos9+Je'<'+Je- 



(3^ 



P 



^{dj - dj) + {Jdj - Jdj) 



Z{J, J) = 

Z{J, J) = 
djdjZ{J, J) = Z{J, J) ; 



dp-^{dj + dj) 



(39) 

(40) 
(41) 



where (39) and (40) are the Schwinger-Dyson and action principle equations, while (41) is a 
constraint equation. In fact these equations have a two parameter (one if you neglect the normal- 
ization) class of solutions given by linear combinations of basis solutions 



Z( J, J) 



^^g/3cos0+Je'^+Je-''' 



(42) 



for contours S equivalent to either Si and S2 in figure 9 (assuming real positive /5). 



— TT 



^1 



TT 









— TT 






^2 







Figure 9. 



Note that integration over Si — S2 is equivalent to the integral over the usual compact path 
6 = [—71, it]. The possible physical relevance of the exotic solutions, upon generalizing to a theory 
in a finite number of dimensions, is not manifest as it was for polynomial theories. In the former 
case, symmetry breaking vacua were clearly related to exotic solutions. Our experience with the 
polynomial theories leads us to speculate that the exotic solutions for lattice gauge theories are 
related to physically realizable phases of gauge theory. Certainly we expect that the appearance 
of phase boundaries in gauge theories, via the accumulation of Lee- Yang zeros, is closely related 
to the collapse of the solution set in the thermodynamic limit. 
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6 Conclusions and Outlook 



We have examined the complete set of solutions of the differential equations which follow from 
the Schwinger action principle. While only one of these solutions corresponds to the usual path 
integral, the other solutions, which involve complexified path integrals have physical relevance. 
On the one hand the manner in which the full solution set collapses in the thermodynamic (or 
analogous qt 0) limit is related, via Stokes phenomena, to the accumulation of Lee- Yang 
zeroes at phase boundaries. On the other hand, exotic solutions may themselves be physical, with 
theta vacua and symmetry breaking vacua as known examples. In the zero dimensional case, we 
have proven that Borel resummations of perturbation series, having various singularity avoiding 
contours of integration in the complex Borel variable, solve the action principle equations and 
therefore correspond to various complexified "path" integrals. 

While we have explicitly discussed the solution set of the Schwinger Dyson and action princi- 
ple equations for zero dimensional models, one can readily generalize to lattice models in multi- 
dimensions, in which case one finds a huge solution set. The basis set of solutions to the Schwinger- 
Dyson equations for a scalar field theory on a lattice can be written in terms of the zero dimensional 
solutions Z^^\j) as follows; 

Z = ^MK.,^^^)\{zf\j,) (43) 

where Kij is the lattice kinetic term, and the zero dimensional solution Z^'^ may be different 
at each lattice site k. For a polynomial potential, the number of independent solutions grows 
exponentially with the number of lattice sites. Like internal symmetries, space-time symmetries 
may be broken by the choice of integration paths. Determining the collapse of the solution set in 
the thermodynamic and continuum limits is a difficult problem. It would be very interesting if 
exotic solutions with different integration paths at different sites have physical relevance. 
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